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We report on the measurement of Stark shifted energy levels of 8 ' Rb Rydberg atoms in static 
electric fields by means of electromagnetically induced transparency (EIT). Electric field strengths 
of up to 500 V/cm, ranging beyond the classical ionisation threshold, were applied using electrodes 
inside a glass cell with rubidium vapour. Stark maps for principal quantum numbers n — 35 and 
n = 70 have been obtained with high signal-to-noise ratio for comparison with results from ab initio 
calculations following the method described in [M. L. Zimmerman et al., Phys. Rev. A 20, 2251 
(1979)], which was originally only verified for states around n = 15. We also calculate the dipole 
matrix elements between low-lying states and Stark shifted Rydberg states to give a theoretical 
estimate of the relative strength of the EIT signal. The present work significantly extends the 
experimental verification of this numerical method in the range of both high principal quantum 
numbers and high electric fields with an accuracy of up to 2 MHz. 

PACS numbers: 32.30.-r, 32.60.+i, 32.80.Ee 


I. INTRODUCTION 

The response of atoms to static electric fields (DC 
Stark effect) results in line shifts, state mixing and, for 
sufficiently large fields, ionisation. The line shifts are con¬ 
ventionally summarised in Stark maps, displaying the en¬ 
ergy levels as a function of the applied field. Stark maps 
of alkali atoms are routinely calculated by diagonalising 
the perturbed Hamiltonian [T] , taking into account quan¬ 
tum defects and corresponding electronic wavefunctions 
[2]. A precise knowledge and control of Stark shifted Ry¬ 
dberg states is required for the application of Rydberg 
atoms as quantum probes [3], for controlling the inter¬ 
actions between Rydberg atoms 5], the production of 
circular Rydberg atoms [5], the structure and dynam¬ 
ics of Rydberg gases mm, and possible applications in 
quantum information processing [Sj. In the context of 
hybrid quantum systems based on atoms and solid state 
quantum circuits mm, Stark shifts of Rydberg states 
EM and their control m are of particular interest. 

Stark maps of Rydberg excited alkali atoms were stud¬ 
ied in the 1970s using pulsed laser excitations and subse¬ 
quent pulsed-field ionisation SHE]- For low-lying Ryd¬ 
berg states of Sodium, Lithium and Caesium with prin¬ 
cipal quantum numbers n < 20, Stark maps have been 
recorded up to and beyond the classical ionisation limit 
[3 US]. Stark maps of 85 Rb for n up to 55 were studied in 
the 1980s for low electric fields using two-photon laser ex¬ 
citation and detecting ionisation from thermal collisions 
m ■ In this regime of low electric fields the Stark maps 


* jens.grimmel@uni-tuebingen.de 
t fortagh@uni-tuebingen.de 


do not show level crossings but quadratic dependence on 
the applied field with slight deviation from this for highly 
excited states (n = 55). 

Many recent experiments on Stark shifts use electro¬ 
magnetically induced transparency (EIT) [T7, IS]. This 
spectroscopic method provides a high resolution of the 
energy levels m and is suitable for the detection of states 
of high principal quantum numbers. However, the mea¬ 
surements so far only covered the range of low electric 
fields, in which just the first avoided crossings appear. 
Similarly, the detection of ions by micro-channel plates 
(MCP) [19] and the technique of measuring the ionisa¬ 
tion currents from Stark-shifted Rydberg states [50] have 
only been used at low electric fields. 

In this article, we report on the optical spectroscopy 
of Stark shifted Rydberg states with principal quantum 
numbers of n = 35 and n = 70 for electrostatic fields 
between 0-500 V/cm, ranging beyond the classical ioni¬ 
sation limit. Our measurements go to higher principal 
quantum numbers and cover a three to four times larger 
range of electric fields relative to the classical ionisation 
threshold than any of the aforementioned works using 
EIT. The observed Stark maps are compared with re¬ 
sults from numerical calculations following the lines of 
the numerical method by [T], including the recently im¬ 
proved accuracy of the quantum defects m- In the orig¬ 
inal work this method was optimised for Rydberg states 
around n = 15 due to computational limitations as well 
as the accessible experimental data. The high accuracy 
of this method at low electric fields is always given by the 
accuracy of the unperturbed energy levels, which are used 
in the calculation. However, for strong electric fields and 
high principal quantum numbers, where a high accuracy 
of the calculated energy levels is desirable, this method 
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FIG. 1. (Color online) Vapour cell with electrodes. The 
capacitor plates are formed by two glass plates with 5nm of 
Nickel facing each other at a separation of 5 mm. Two wires 
are glued to each capacitor plate to contact them. For the 
measurements the laser beams counterpropagate through the 
view ports indicated in the image. 


becomes more susceptible to numerical errors jT]. In this 
article we show by direct comparison, that it is still ap¬ 
plicable even in these regions of study. Furthermore, we 
calculate dipole matrix elements between 5P 3 / 2 and the 
observed Stark shifted states, which are then used to give 
an estimate for the relative strength of the measured sig¬ 
nals. 


II. MEASUREMENT OF STARK MAPS 


For the measurement of Stark shifts of 87 Rb Rydberg 
atoms we use a vapour cell with a pair of plate electrodes 
for applying homogeneous electric fields. The electrodes 
are inside the cell [Fig. [l], trying to avoid the effect ob¬ 
served in previous works where static electric fields ap¬ 
plied to a vapour cell with outside electrodes are com¬ 
pensated by ionised Rubidium and electrons [22]. The 
Rubidium vapour in the cell is at room temperature with 
an estimated pressure of ~ 1 • 10~' mbar. The electrodes 
in the center of the glass cell [Fig. [I] are formed by two 
square glass plates coated with 5 nm of Nickel. This gives 
the plates ss 60 % transparency, which allows for optical 
access on the axis perpendicular to them. However, this 
optical access was not used in the experiments presented 
in this work. The plates are mounted with 5 mm sepa¬ 
ration on insulating ceramics (Macor) and connected to 
a voltage source through metallic wires. Laser beams 
for the optical spectroscopy are introduced through view 
ports and pass through the cell between the plates. 

We measure the transition frequency from the Ru¬ 
bidium ground state to Rydberg states by electromag- 
netically induced transparency (EIT). The EIT three 
level ladder scheme [22] consists of the ground state 
5Si/2(T’=2), the intermediate state 5P 3 / 2 (F=3) and a 
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FIG. 2. (Color online) Ladder scheme for electromagneti- 
cally induced transparency (EIT). For our measurements the 
probe laser is held on resonance with the 5 Si/ 2 (-F’=2) —> 
5 P 3 / 2 (E=3) transition of s 'Rb and the coupling laser is 
scanned around the upper transition. Dipole selection rules 
only allow for the measurement of nS and nD states in this 
scheme. 


Rydberg state nS or nD [Fig. [ 2 ]. However, it is im¬ 
portant to keep in mind that l is not a good quantum 
number any more in the presence of an external elec¬ 
tric field and that the Stark shifted states can be con¬ 
sidered as a mix of all possible unperturbed l states. 
For the spectroscopy we use a probe laser that is locked 
to the 5Si/ 2 (F=2) -a 5P 3 / 2 (F=3) transition (780nm) 
and a coupling laser with a variable frequency close to 
the transition between 5P 3 / 2 (F=3) and a Rydberg state 
(480nm). When the coupling laser is resonant with this 
transition we detect a maximum in the transmission of 
the probe laser through the vapour. 

We use a grating-stabilised diode laser (Toptica, DL 
pro) of ss 100 kHz linewidth as the probe laser and 
a frequency-doubled, grating-stabilised diode laser of 
ss 200 kHz linewidth (Toptica, TA-SHG pro) as the cou¬ 
pling laser. Both lasers are phase-locked to a frequency 
comb (Menlo Systems, FC 1500). For conveniently se¬ 
lecting the right modes of the frequency comb for both 
lasers we use a calibrated wavelength meter (HighFi- 
nesse, WS Ultimate 2) [21]. The power of the probe laser 
is 1 pW and the power of the coupling laser is 25 mW 
with 1/e 2 diameters of 450 pm and 150 pm in the cell, 
respectively. The small diameter of the coupling beam 
results in a high intensity and therefore high Rabi fre¬ 
quency on the corresponding transition while the bigger 
diameter of the probe beam is chosen in order to ensure 
maximal overlap of the two beams within the cell. The 
laser powers and polarisations were adjusted to maximise 
the EIT signal in zero field [Fig. [ 3 ]. A measurement of 
the polarisations after this adjustment revealed that this 
resulted in circular polarisations for both lasers. The 
frequency of the probe laser is kept on resonance with 
the 5Si/ 2 (F=2) -A 5P 3 / 2 (E=3) transition [Fig. [ 2 ] and 
its transmission through the cell is measured using an 
avalanche photodiode. 

In order to reach a high signal-to-noise ratio in the 
spectroscopy, we use the following lock-in measurement 
technique. We modulate the power of the coupling laser 
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FIG. 3. (Color online) Optical setup (schematic) for the mea¬ 
surement of Stark maps using EIT. Both lasers are guided to 
the experiment through optical fibres and overlapped within 
the volume of the plate capacitor using two dichroic mirrors. 
Intensity and polarisation of both lasers are adjusted inde¬ 
pendently. The transmission of the probe beam through the 
vapour cell is measured with a photodiode. 


using an acousto-optic modulator (AOM) and demodu¬ 
late the signal from the photodiode on the carrier fre¬ 
quency. In addition we modulate the frequency of the in¬ 
tensity modulation and demodulate the photodiode sig¬ 
nal on two of the sidebands that arise from this. We 
then average the demodulated signal from the carrier 
and both sidebands. The additional frequency modu¬ 
lation decreases the signal strength on the carrier fre¬ 
quency and creates an even lower signal strength on the 
sidebands, but averaging the signal from three frequen¬ 
cies lowers the noise level at the same time. In total, the 
resulting signal-to-noise ratio is improved by a factor of 
~ 2 as compared to a simple lock-in measurement with¬ 
out frequency modulation. For every set frequency of the 
coupling laser we ramp up the voltage on the plate ca¬ 
pacitor using an auxiliary output of the lock-in amplifier 
(Zurich Instruments, HF2LI). This is sufficient for elec¬ 
tric fields up to ss 20V/cm. For even higher fields up to 
~ 500 V/cm we use an additional voltage amplifier. 

Measured data near the unperturbed 35S;iy 2 state is 
shown in Figs. [4] and [5j Near the unperturbed 70S!/ 2 
state we conducted a preliminary measurement, which is 
not shown in this work, and then selected a smaller re¬ 
gion for a more detailed measurement in order to test the 
frequency precision of the numerical calculations [Fig. [6]. 
Each pixel in the gray scale images represents an average 
of the demodulated signal over 200 ms. Between every 
two pixels we add a waiting time of 50 ms to allow for 
the low-pass filter of the lock-in amplifier to settle. The 
gray scale was adjusted with a cutoff for better visibility 
of weaker signals. Our data shows states ranging up to 
and even beyond the classical ionisation threshold that 
have not been measured by means of EIT before. The 


classical ionisation threshold Ej on , i.e. the saddle point 
which is formed by a Coulomb potential with an external 
electric field F, is given by 

E ion = -2 VF (1) 

in atomic units [2]. This results in an electric field 
strength for the ionisation threshold of ss 312 V/cm for 
35 Si/ 2 and ~ 16 V/cm for 70 Si/ 2 . For a quantitative 
analysis we give a brief review of the numerical calcula¬ 
tion of Stark maps in section |III| 

We observe two background effects [Fig. [4j [6] which 
are caused by the region of the cell that is not covered by 
the plate capacitor [Fig. [lj. The first is the line of the 
unperturbed state, which remains visible for all applied 
voltages because we probe those outer regions of the cell 
as well. The second is a smearing of the lines to the right 
at low fields as visible in Fig. [4] a t the avoided crossings 
up to ss 50 V/cm and in Fig.d6l For a certain electric 
field strength inside the plate capacitor one always finds 
lower electric field strengths in the inhomogeneous out¬ 
side region, causing the asymmetry of the smearing to 
the right. 

Inside the capacitor undesirable electric fields could 
also arise from the dipole which is formed between ad¬ 
sorbed Rubidium and the Nickel surface of the capacitor 
plates. In previous experiments the repeated deposition 
of cold atom clouds of Rubidium onto a Copper surface 
led to electric fields close to the surface which saturated 
as the number of deposited clouds increased HU. Since 
the work functions of Nickel and Copper are similar and 
the cell is filled with Rubidium vapour at all times we 
may find a similar effect for the capacitor plates. This 
effect may play a role for Rydberg states higher than 
n = 70 at low electric fields, but is negligible for the 
measured data presented in this article, where the dis¬ 
tance between the laser beams and the capacitor plates 
is « 2 mm. Other stray fields could originate from out¬ 
side of the cell, but should be compensated by electrons 
and ionised Rubidium [22] . 


III. CALCULATION OF STARK MAPS AND 
DIPOLE MATRIX ELEMENTS 

The Hamiltonian for the valence electron in an alkali 
atom with an external electric field in z-direction can be 
written as 

H = H 0 +Ez (2) 

in atomic units with Hq the Hamiltonian for the valence 
electron in absence of any perturbation, E the electric 
field strength and z the position operator in z-direction. 
For the following calculations we include the fine struc¬ 
ture splitting in Ho and neglect the hyperfine structure 
splitting as it is smaller than our frequency resolution in 
the experiment [T] . The Stark shifts of the 5S!/ 2 (F=2) 
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FIG. 4. (Color online) (a) Stark map for 35 Si/2 for an electric field range of 0-235 V/cm. The gray scale is determined from 
the measured EIT signal, which represents the demodulated transmission signal from the photodiode. Due to the frequency 
resolution of 125 MHz some of the lines with a small slope appear dotted here, (b) The red, green and blue lines show the 
numerically calculated Stark map for \rrij\ = 1/2, \rrij\ = 3/2 and \rrij\ = 5/2 respectively. The opacity gradients indicate 
the calculated transition strength D to the corresponding Stark shifted Rydberg states. In this region we find a very good 
agreement between experimental and numerically calculated results in the Stark shifts as well as the transition strengths. 


and 5P 3 / 2 (.F=3) states from the EIT scheme are negligi¬ 
ble as well, as they amount to less than 100 kHz for the 
region of the field strength up to 500 V/cm. We then cre¬ 
ate a matrix representation of H in a subset of the basis 
given by Hq. This way Hq is represented by a diagonal 
matrix with the energy levels from m on its diagonal. 

The matrix representation of z is symmetric with only 
off-diagonal entries and in spherical coordinates we ob¬ 
tain from pQ 


{n,l,j,mj | z | 

= i {n,l,j\r\ri,l',j’) 


Y | \' mu ~ mi 


mi — rn j =t k 


1 


X 2 ,rni,mj ~ m 
x (l, mi | cos 0 1 1', mi) 


J ,rn 3 


(3) 


with the radial overlap integral (n, l, j \ r \ n !, l',j') in the 
first line of the right hand side. From the Kronecker delta 
we can see that the matrix representations of H 
for different values of m 3 can be calculated separately. It 
is also evident that the Stark shifts for ±|wj| are always 
equal. Therefore we only calculate matrix representa¬ 
tions for positive values of rn 3 in this step. This reduces 
the computing time for diagonalising the matrices, but 
we need to consider states with both signs for the calcu¬ 
lation of transition strengths later on again. 

For the calculation of the overlap integrals we tried 
two different methods to calculate the radial wavefunc- 
tions for all states in the chosen subset of the basis and 
also for 5 Si/ 2 (.F=2) and 5P 3 / 2 (F=3), which will be used 
later on in the calculation of the measure for the tran¬ 
sition strengths. One method was to further follow |I] 
using their atomic potential for 87 Rb and Numerov’s 
method to solve the radial Schrodinger equation. We 
also tried using the parametric model potential from [23] 
with another solving algorithm [23], which allowed us to 
obtain parts of the wavefunctions that are located fur- 
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FIG. 5. (Color online) (a) Stark map for 35 Si/2 for an electric field range of 235-470 V/cm. The gray scale represents the 
measured EIT signal. Due to the frequency resolution of 125 MHz some of the lines with a small slope appear dotted here. The 
blue solid line indicates the classical ionisation limit. Beyond this limit fewer states are present in the measurement. The states 
also appear broader and weaker as the external electric field lowers and extends the opening of the potential barrier it forms 
with the atomic potential, (b) Results from our numerical calculations for the same region. The red lines represent \rrij\ = 1/2, 
the green lines \m,j\ =3/2 and the blue lines \m.j\ =5/2 with the opacity gradients indicating the calculated transition strengths 
D to each state, (c) and (d) show details of the measurement with the numerically calculated Stark map on top. The opacity 
of the colors has been increased with a linear scaling in comparison to (b) for better visibility of the calculated lines on top of 
the experimental data. The agreement in the energy levels remains even beyond the classical ionisation threshold, but we find 
some discrepancies in the transition strengths as the electric field strength increases. 
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FIG. 6. (Color online) Detailed part of a Stark map for 
70 Si/ 2 for an electric field range of 6.23-6.70 V/cm. The gray 
scale represents the measured EIT signal. The red, green and 
blue lines plotted on top of the measured signal show the nu¬ 
merically calculated Stark map for \rrij\ = 1/2, \rrij\ = 3/2 
and \rrij\ =5/2 respectively with the opacity gradient scaled 
to the calculated transition strength D. The frequency resolu¬ 
tion of this measurement is 2 MHz and, again, we find a very 
good agreement between the experimental and numerically 
calculated results. 

ther inside the ionic core. While the two methods lead 
to slightly different amplitudes of the wavefunctions due 
to normalisation, the differences in the Stark shifts and 
transition strengths calculated from the two methods for 
the observed states lie below our experimental accuracy. 
All calculations presented in this work utilised the latter 
method. 

We can then calculate the matrix representations of 
the perturbed Hamiltonians H using equation ([3| and 
diagonalise them efficiently in parallel for different values 
of the electric field strength to obtain the energy eigen¬ 
values. In this step we also calculate the eigenvector /3 
corresponding to every eigenvalue. With this we can fur¬ 
ther follow the method from p] to represent every Stark 
shifted Rydberg state £ as a linear combination of unper¬ 
turbed states, given by the eigenvectors, and calculate 
the dipole matrix elements 

(£ | r | 5P 3/2 , F=3, mp) 

= 53 Pn,i,j,\m j \(nJ,j,m j \r\5P 3 / 2 ,F=3,m F ) 

n : l,j,rrij 

( 4 ) 

with /3 nt ij t \ m .\ the entry of /3 that corresponds to the 
states \n,l,j,i\nrij\). Here we sum over all states in the 
selected subset of the basis given by the unperturbed 
states and both signs for each value of rrij. For the cal¬ 
culation of the dipole matrix elements between the un¬ 
perturbed Rydberg states and |5P 3 / 2 , F=3, m F ) we sep¬ 
arate the radial and the angular parts of the wavefunc¬ 
tions. We use the radial overlap integrals from equation 
[3] and calculate the three components of the angular over¬ 
lap integrals following [251 ]. 


We take into account the effect of the probe laser, 
which couples the different m? substates of |5 Si/ 2 , F=2) 
and |5P 3 /2,F'=3) by calculating the weighting factors 

Vm p = 5Z I ( 5P 3/2,-F=3, m F | £ p f I 5Siy 2 , F'=2, m' F ) |~ 

m' F 

( 5 ) 

with £ p the polarisation of the probe laser. In this step we 
assume that all mj? substates of |5 Si/ 2 , U=2) are evenly 
occupied. The weighting factors are then used to calcu¬ 
late the measure for the transition strength 

D = 5^ Vm F I (£ I £ c r\ 5P 3 /2, F= 3, m F ) \ 2 (6) 

m f 

with e c the polarisation of the coupling laser. In the 
experiment the two lasers are counterpropagating and 
perpendicular to the external electric field with circular 
polarisations. Therefore we use £ p = e c = (0,1,1) here. 
The measure for the transition strength D is used to de¬ 
termine the opacity gradients for the lines in Fig. 00 
and [b] 

The Stark map for 35 Si/ 2, which is shown in Figs. [4] 
and [5j was calculated from a subset of the basis of ap¬ 
proximately 1600 states and for 2000 values of the electric 
field strength using MATLAB. For 70S!/ 2 ) shown in Fig. 
[6] the calculations run similarly with 4000 states and 500 
values of the electric field strength. The number of states 
that was used for these calculations was adjusted so that 
a further increase only yields changes which lie below 
the accuracy of the experimental data. Calculated Stark 
maps for \rrij\ = 1/2, \rrij\ = 3/2 and \rrij\ = 5/2 are in¬ 
cluded in the figures. The measurements were performed 
in frequency regions around S states, but since l is not 
a good quantum number any more in the presence of an 
electric field pQ, we find that other states with rrij \ = 3/2 
and rrij = 5/2 are shifted far enough by the Stark effect 
to appear within the measured frequency range. 


IV. COMPARISON OF MEASUREMENTS AND 
CALCULATIONS 

Figs. [4] and [5] show comparisons of the measured and 
numerically calculated Stark maps for 35Sx/2- The fre¬ 
quency axis shows the detuning of the coupling laser rel¬ 
ative to the absolute value of the transition frequency 
from [2T|. A linear scaling with an offset has been ap¬ 
plied to the electric field axis of the measured data. Using 
only this scaling we achieve a match between calculated 
and measured energy levels for Stark states in the whole 
range of our measurement. 

The lines in the calculated Stark maps are drawn with 
an opacity gradient, resulting in a color range between 
white and the respective color associated with the dif¬ 
ferent values for rrij \. The opacity gradient is scaled 
to D from equation ([6]) with an upper cutoff at 70% 
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of its maximum for better visibility. We find a good 
agreement between this calculated measure for the tran¬ 
sition strength and the experimental data in the range 
of electric fields shown in Fig. [4j The quality of the 
agreement deteriorates slightly approaching the classical 
ionisation threshold, but there is still good agreement 
beyond this point, as can be seen in Fig. [5j The very 
weak signals we still find in the experimental data at 
electric fields around 400V/cm all correspond to calcu¬ 
lated states with a strong calculated transition strength. 
Considering this, we think that the method described 
here can be applied for these high electric field strengths, 
where only few states are still visible in the experimen¬ 
tal data, as a way to select potentially interesting areas 
before a measurement. 

One important cause for discrepancies between the ex¬ 
perimental and numerically calculated data we present 
here is that no ionisation effects were taken into account 
for the calculations. However, it is interesting to note 
that even though the calculated measure for the tran¬ 
sition strength shows some differences at high electric 
fields, we still find a remarkable agreement of calculated 
and measured energy levels in this region of study. These 
differences could be related to the calculation of the ra¬ 
dial wavefunctions. We calculate the radial wavefunc- 
tions for |5S]y 2 ) and |5P 3 / 2 ) using the same method as 
for Rydberg states. Evidently this works quite well here, 
but more accurate radial wavefunctions, especially for 
|5P 3 / 2 ) but also for Rydberg states with low l, could 
further improve the quality of the calculated transition 
strengths. This is further emphasised as the biggest dis¬ 
crepancies of the transition strengths at high electric 
fields can be found for states \rrij\ = 1/2 while states 
with \rrij\ = 5/2 show hardly any discrepancies on the 
whole range of electric fields. On the other hand, these 
discrepancies could also be explained by a stronger cou¬ 
pling of Rydberg states with low l to the continuum. 
Some disagreements on the whole range of our measure¬ 
ments can also be caused by the scaling of the EIT signal 
with the dipole matrix elements, which is only approxi¬ 
mated by the squared scaling in D. Furthermore, lines 
with small slopes and especially local extrema may ap¬ 
pear over-pronounced in the experimental data. Since 
these parts of the lines are broad in the electric field do¬ 
main and the inhomogeneous electric field outside of the 
plate capacitor still lies within the volume of the cell, this 
results in an effective contribution of more atoms to the 
signal. 

We find a match similar to the one presented for 35S]y 2 
in the energy levels and transition strengths for 70S;i/ 2 , 
which is shown in detail in Fig. [G] Discrepancies lie 
within the frequency resolution of 2 MHz for Fig. [6j 
which is the highest accuracy we present here. The opac¬ 


ity gradients were assigned the same way as for 35 Si/ 2 . 
The calculated transition strengths match the experi¬ 
mental data equally well as for 35S;iy 2 . The experimental 
data presented in Figs. 00 and [6] is provided as supple¬ 
mentary data to this article [2T>/ All results presented 
in this work were measured and calculated in the energy 
range near nS states. However, they intrinsically include 
contributions all possible l states due to the mixing of 
states caused by the external electric held. 

For a possible extension of our measurements to a 
wider range of principal quantum numbers n than the re¬ 
gion between n = 35 and n = 70, the most important lim¬ 
iting factor is that the coupling of Rydberg states to the 
intermediate 5P 3 / 2 (.F=3) grows weaker as n increases, 
which leads to a lower signal-to-noise ratio. Another fac¬ 
tor influencing the signal-to-noise ratio is the density of 
the vapour, which on the other hand gives rise to collec¬ 
tive effects which in turn cause undesired changes of the 
EIT signal [27]. 

V. CONCLUSION 

In summary, we have demonstrated agreement between 
measured and calculated Stark maps of ®'Rb up to an 
accuracy of 2 MHz. This agreement holds for the range 
from zero held to beyond the classical ionisation thresh¬ 
old and for principal quantum numbers n = 35 and 
n = 70. Such data and calculations may aid the accurate 
mapping of electric helds at surfaces m- Furthermore, 
we presented numerical calculations to estimate the tran¬ 
sition strength from low-lying states to Stark shifted Ry¬ 
dberg states. The results from these calculations show 
a very good agreement with our experimental data on a 
high range of electric helds and even beyond the classical 
ionisation threshold. 

Altogether, the methods presented in this paper can be 
used to hnd experimentally accessible Stark shifted Ry¬ 
dberg states with an appropriate sensitivity to external 
electric helds for a wide range of applications. For exam¬ 
ple, the dependence of the transition strength to different 
Stark shifted Rydberg levels could be used to determine 
not only the strength of an external electric held in an 
experimental system but also the electric held axis. 
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